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摘 E: 针对 在 室内 环境 的 轨迹 跟踪 过 程 中 ， 通 信 开 销 较 大 和 算法 复杂 度 较 高 的 问题 ， 对 一 种 基于 CSI(channel state 
information) 信 号 的 人 员 轨 迹 跟 踪 方法 进行 研究 。 首 先 从 CSI 中 提取 代表 目标 位 置 (角度 ) 概 率 的 AOA (angle-of-arrival) 9t 
谱 ， 通 过 MUSIC 算法 得 到 的 多 普 勒 频 移 与 AOA 频谱 相 结合 来 确定 人 员 移 动 速度 及 位 置 ; 最 后 利用 改进 的 三 边 定位 质 
心算 法 确定 人 员 位 置 ， 模 拟人 员 移 动 轨迹 ， 实 现 对 室内 人 员 精 准 的 跟踪 定位 。 通 过 与 其 他 算法 以 及 不 同人 员 移 动 速度 
的 对 比 ， 仿 真实 验 表明 本 文 提出 的 人 员 追 踪 方法 可 以 大 幅度 提高 定位 的 精确 性 与 稳定 性 
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Abstract: For high communication overhead and greater algorithm complexity of the problem, researchers on a progressive 


track CSI (channel state information) signal tracking method based on the trajectory of the indoor environment of the tracking 


process. Extracting first representative of the CSI target position (angle) probability AOA (angle-of-arrival) spectrum obtained 


by the Doppler shift Music AOA algorithm in combination with spectral determine the moving speed and position the art; finally 
the improved trilateral positioning centroid algorithm to locate personnel, modellers moving track to achieve accurate indoor 
tracking and locating personnel. By comparison with other algorithms and different speeds of movement of persons, person 


tracking simulation results show that the proposed method can greatly improve the accuracy and stability of positioning. 
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没有 考虑 到 多 人 追踪 、 实 时 跟踪 、 多 径 效 应 的 影响 问题 。 文 献 

[2] 使 用 软件 无 线 电 (USRP) 实 现 了 整个 实验 区 域 的 手势 识别 ， 

定位 技术 作为 普 适 计算 和 物 联网 的 一 个 重要 组 成 部 分 ， 受 使 用 FFT 来 检测 多 普 勒 频 移 。 文 献 [3] 中 使 用 RFID 构建 虚拟 触 

到 了 越 来 越 广泛 的 关注 。 全 球 卫星 导航 系统 (global navigation R, 并 利用 AoA 信息 进行 细 粒 度 跟踪 。WiDraw 利用 AoA 信 

satellite system，GNSS) 在 室外 精准 定位 方面 起 着 关键 作用 , 但 息 很 好 的 提供 了 使 用 相位 信息 确定 位 置 的 思路 ， 但 需要 Wi-Fi 

是 由 于 在 室内 环境 下 ，GNSS 由 于 信号 在 传播 过 程 中 表现 出 来 ”发 射 机 在 环境 中 覆盖 所 有 方向 内 。 而 在 实际 环境 中 受 多 径 效 应 

的 多 径 效 应 、 散 射 、 衍 射 等 特性 ， 使 其 无 法 再 发 挥 其 高 精度 定 ”的 影响 ， 实 现 起 来 存在 一 定 难度 。 文 献 [5] 中 利用 CSI， 构 建 了 

I 优势。 目前 基于 CS 的 定位 技术 研究 已 经 取得 了 一 定 的 成 基于 CSI 模型 的 低 成 本 高 精度 被 动 式 目标 定位 方法 LIFS， 将 
果 ， 甚 至 已 经 达到 了 米 级 或 是 亚 米 级 的 定位 精度 。 CSI. 的 特点 有 效 的 结合 在 了 目标 定位 方面 ， 但 未 考虑 检测 区 二 

近 几 年 以 来 ， 麻 省 理工 学 院 、 斯 坦 福 大 学 、 华 盛 顿 大 学 、 与 检测 率 的 关系 ,文献 [6] 中 将 CST 中 的 相位 信息 应 用 在 通过 人 

香港 科技 大 学 、 清 华 大 学 等 高 校 都 对 CSI 感知 应 用 系统 进行 了 员 说 话 时 口腔 的 运动 ， 很 好 的 达到 更 细 粒 度 的 定位 的 效果 ， 但 

大 量 的 研究 。 文 献 [1] 中 使 用 商用 Wi-Fi 实现 无 设备 室内 人 员 跟 ” 是 相位 信息 不 易 获 取 ， 并 且 要 利用 专用 设备 USPR， 稳 定性 和 

踪 ， 实 现在 丢 包 、 延 迟 的 情况 下 也 能 实现 很 高 的 定位 精度 ， 但 “实用 性 不 高 。 文 献 [7] 中 FIMD 系统 通过 利用 CSI 的 稳定 性 ,在 
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态 环境 下 实现 了 更 细 粒 度 的 人 员 检 测 ， 但 未 达到 高 检测 率 ， 
并 且 会 受 实验 环境 的 改变 而 影响 整体 性 能 ,文献 [8] 中 的 BFP 系 
统 利 用 CSI 进行 行为 无 关 移 动 检测 ， 整 体 性 能 较 好 ， 但 其 算法 
的 整体 效率 不 高 。 

针对 现 有 的 轨迹 跟踪 方法 大 都 是 以 牺牲 较 高 的 通信 开销 和 
较 高 算法 复杂 度 为 代价 来 提高 精度 ， 本 文 提 出 的 基于 CSI 信号 
的 人 员 轨 迹 追 踪 方 法 很 好 地 规避 了 上 述 存在 的 问题 。 首 先 提取 
到 的 CSI 信号 用 平滑 算法 进行 处 理 ， 形 成 增强 型 的 CSI 算法 。 
增强 型 CSI pun 前 一 个 CSI 的 值 大 小 与 当前 的 CST 值 进行 
对 比 ， 降 低 了 移动 人 员 目 标的 动能 波动 。 最 后 再 利用 改进 的 三 
边 定 位 质心 算法 确定 每 一 个 活动 点 ， 将 所 有 的 活动 点 绘制 成 轨 
迹 图 像 直观 显示 出 来 ， 实 现 对 室内 人 员 进 行 更 加 精准 的 跟踪 定 


位 。 
1 ”相关 理论 


1.1 信道 状态 信息 
CSI 是 衡量 信道 情况 的 信道 状态 信息 ,属于 PHY 层 ， 来 自 
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bE, 
1.2 CSI 信号 的 多 普 勒 频 移 提取 
多 径 传播 就 是 在 带 有 一 对 发 射 器 和 接收 器 的 室内 环境 中 ， 
言 号 不 仅 沿 直达 路 径 传播 ， 还 存在 被 其 他 物体 (如 人 ) 和 墙壁 让 
挡 反 射 的 现象 0553。 在 接收 器 处 接收 的 信号 因此 是 来 自 所 有 路 
径 的 信号 的 做 加 。 当 一 个 人 在 环境 中 移动 时 , 路 径 的 长 度 受 人 
的 反射 信号 不 同 相 应 地 改变 ， 这 时 引起 载波 频率 多 普 勒 频 移 反 
射 信号 发 生 改变 ， 载 波 频 率 如 式 (6) 所 示 。 


f Doppler 


= faet (6) 
其 中 : 了 是 信号 的 原始 载波 频率 ,vw 是 路 径 长 度 变化 的 速度 ， 
c 是 Wi-Fi 信号 在 空中 的 传播 速度 。 当 人 在 移动 时 ， 引 入 多 普 
P Wi-Fi 信号 的 频率 转换 对 于 5GHz ERUR H JLI. 
显然 对 细 粒 度 如 此 小 的 多 普 勒 载波 频 移 进行 检测 是 十 分 困难 的 。 
2 基于 CSI 的 轨迹 追踪 方法 


2.1 基于 MUSIC 的 多 普 勒 估计 
状态 信息 在 采集 时 ， 由 于 实际 环境 噪声 和 干扰 造成 的 分 组 


OFDM 系统 下 解码 的 子 载波 [9-10]。CSI 是 细 粒 度 的 物理 信息 ， 
对 环境 更 加 敏感 ， 应 用 于 动作 识别 ， 手 势 识 别 ， 击 键 识别 ， 跟 
踪 等 领域 。 首 先 需要 知道 被 测 节点 的 收 功率 ， 发 射 节 点 信号 功 
率 作为 参考 ,然后 利用 信号 传播 衰减 模型 转换 为 节点 间 的 距离 ， 
最 后 再 用 改进 的 三 边 测量 算法 就 可 以 简单 确定 未 知 节点 的 位 置 
无 线 信 号 的 发 射 功率 和 接收 功率 之 间 的 关系 可 以 ) 
P, - P, 1d" a) 

其 中 : 及 是 无 线 信号 的 接收 功率 , P 是 无 线 信号 的 发 射 功率 , 
是 收发 单元 之 间 的 距离 , 传播 因子 ,数值 大 小 取决 于 无 线 信号 
传播 的 环境 。 
对 式 (1) 的 两 边 取 对 数 可 以 得 到 式 (2)。 


10*nlg r - 10g Q) 


节点 发 射 功率 是 已 知 的 ,将 发 送 功率 代入 式 (2) 中 得 式 (3)。 
l0lg P, = A-10*nlgr 6) 
式 (3) 的 左 半 部 分 10lg 玉 是 接收 信号 功率 转换 为 4Bm 的 表 
达 式 ,可 以 直接 写成 式 (4), 在 式 (4) 中 A 可 以 看 做 信号 传输 lm 远 
时 接收 信号 的 功率 . 
P.(dBm) = A-10*nlgr (4) 
式 (4) 中 有 (dBm) 就 是 CSI 常数 ，A 入 的 数值 决定 了 接 
收 信号 强度 和 信号 传输 距离 的 关系 。 
因为 每 个 参考 点 的 信号 衰减 程度 不 同 ， 所 以 只 有 一 个 1 的 


值 不 具有 一 般 性 ， 也 很 难说 明 其 准确 性 。 由 式 (4) 得 到 式 (5)， 从 
而 得 到 一 组 n 的 值 ， 
NE AS T 
101g 


GR Hi B5] —28 n; ROTE mms ns n, RIX 2H AV 3C 


n» n 这 里 的 即 为 所 求 的 传播 因子 。 


] 式 (1) 表 示 。 


丢失 /延迟 ， 容 易 引 起 不 能 正常 发 送 、 接 收 样本 数据 包 的 问题 。 
为 了 解决 这 些 问 题 ， 如 何 使 用 Wi-Fi 设备 获得 准确 的 多 普 勒 频 
移 估 计 成 为 关键 ， 本 文 提出 了 一 个 基于 MUSIC 的 算法 来 获得 
准确 的 多 普 勒 频率 转移 估计 的 算法 。 
设 在 加 时 刻 收集 第 一 个 CSI 样本 ， 每 个 样本 对 第 一 个 样本 
的 采样 间隔 为 [0,Ab，…Atx] ,其 中 at; 70 。 在 一 个 短 采样 窗 
因为 不 同 CSI 样 JS RDERGESHERL, rbi fne 
Vyan 视 为 一 个 常数 。 第 i 个 CSI 之 间 的 相位 差 样本 ， 第 一 个 样 


号 的 原始 载波 频率 。 因此 ， 阶 段 第 


vArj : a 
Jk: ， 其 中 上 了 是 信 


Mth 个 CSI 样本 与 第 一 个 CSI 样本 之 间 的 差异 可 以 表示 如 下 : 


vA 


a(v) 去 [L, e |f f eI? es e |f F (7) 


其 中 ， go) 多 普 勒 矢量 和 具有 M 个 样本 的 CSI 样本 矩阵 表示 


XCf) -IxCf t). xCf .t; dios 


vA vA 


,© -j2zf * 


Xf ty + Aty )]= 


[L e -j2zf * - Qu 
> 


q x(f,to) +n(f) 


-a(v)x(f.t,) 4n(f) (8) 


Mp: nO 是 噪声 。 当 只 有 一 个 路 径 信号 存在 时 ， 可 以 很 容易 
地 从 CSI 样本 的 相位 测量 中 计算 多 普 勒 频 移 。 在 具有 多 路 径 的 
实际 情况 中 , 工 路径 信号 将 到 达 接 收 器 。 由 (7) 式 和 式 (8)CSI FE 
本 矩阵 可 表示 为 


X(f) =Y atoxt.t,)*n(b = 
i-i 


> > E - AS(f) - NCf) (9) 
[a(v), a(v;).---, a0 Es Cf. f). 


S4 Cf tS). STA S,CÉ S x(t) -N(f) 
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væ i 


A = [a(v,), av, 


普 勒 矢量 是 在 第 一 


个 路 径 信 号 的 路 径 变 化 速度 ; 


a(v,] 是 全 部 的 M* 工 矩阵 多 普 勒 频 移 


个 采样 时 间 测量 的 第 i 个 路 径 信 号 的 CSI, 


噪声 矩阵 。 
个 CSI 样 
载波 上 提供 


T 


x (Ft At) [X(f it +At;) x (E, ,t+AG) X 


FERE SCA) Es Cf 5). so Cf 19)» S, CFT. 是 信号 矩阵 
为 了 得 到 每 个 路 径 信和 号 
本 频率 域 中 的 多 个 快照 ， 如 图 
CSI. 设 K 是 CSI 的 子 载波 数 ， 
CSI 样本 。 对 于 第 i 个 CSI 样本 ， 


, NO 
的 多 普 勒 频 移 。 本 文采 取 每 
1 所 示 。 设 备 在 多 个 子 
有 CSI 的 天 个 快照 


可 以 表示 为 


I 


,tot+Ati)] 00) 


其 中 : 大 是 第 天 个 子 载波 的 原始 载波 频率 。 信 和 号 矩阵 SC/) 中 的 


郝 占 军 ， 等 


XC) SIxCf t) Cf if At), 


M 个 特征 


其 中 : 
方差 。 


对 应 于 噪声 ， 而 其 他 工 特征 值 对 应 于 L 路 径 信号 。 R 


相关 矩阵 Rx 
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xf ty + A 


-[a(v,) av), Ay) 
*[s Cf 15), 5, Cf 15), s, Cf tO) Et NE) 
-AS(f) +N(f) (12) 
] MUSIC 算法 来 估计 使 用 CSI 样本 的 多 普 勒 频 移 。 具 有 
值 的 相关 矩阵 Ry 可 以 表示 为 
Rx = E[XX" ]- AE[SS"]A" + E[INN"] 
=AR A" +o 03) 
s 是 信号 矩阵 的 相关 矩阵，7 是 单位 矩阵 ，5? 是 噪声 的 


具有 M 个 特征 值 。 最 小 的 M -区 特征 值 


第 i 个 路 径 信号 可 表示 为 对 应 于 最 小 的 M-L 特征 值 构 造 一 个 噪声 子 空间 
SU) m Estes rh) S ( fist)] QD — E,-[e...e, ,]。 信 号 和 噪声 子 空间 是 正 交 的 , 所 以 多 普 勒 束 
其 中 : sko) 是 在 第 一 采样 时 间 (4) 第 i 个 路 径 信 号 的 第 个 ” 度 谱 函 数 可 以 表示 为 
快照 的 CSI。 为 了 提高 多 普 勒 频 移 估计 的 准确 性 ， 从 所 有 的 子 NN 1 (14) 
载波 中 选取 CSI， 从 而 避免 选择 低 信 噪 比 的 特定 子 载波 。 由 式 a ()-E,Ey atv) 
(10) 和 (11) , X (9) 可 以 被 表示 为 其 中 : 尖峰 对 应 频谱 函数 的 最 大 值 处 。 由 于 信号 的 原始 载波 频 
率 和 传播 速度 已 知 ， 可 以 获得 多 普 勒 速度 
一 人 K 个 子 载波 ->K 个 快照 
[snis zd CSla CShə . CSl 
一 人 人 一 M 个 数据 包 
CSb; CSLh, . CSLh. 
mui M 个 样本 
pcs CSIm! CSlus — CSImx 
发 送 端 接收 端 Doppler-Music HJ CSDE RE 
图 1 从 M 个 分 组 中 提取 每 个 样本 个 副 载 波 的 KK 个 CST RER 
2.2 ”改进 的 三 边 定位 质心 算法 
ie f PEET , =J a-y (15) 
在 复杂 室内 定位 的 时 候 ， 检 测 到 的 人 员 移 动 轨迹 是 是 时 刻 E 
变化 的 0419， 再 加 上 外 界 不 确定 因素 的 影响 。 导 致 得 到 的 CSI 为 了 避免 侦 然 误差 ， 由 式 (15) 得 到 式 (16)。 
言 号 的 值 波动 较 大 ， 从 而 导致 对 人 员 的 定位 产生 较 大 的 误差 ， 1 
i n d(x,y) 2 -C(dGg, y) -dGo. y) ++ d(x,y, D — (16) 
达 不 到 预期 的 效果 ， 本 文 引入 平滑 算法 ， 把 相 邻 的 CST AHA n 
进行 对 比分 析 ， 去 掉 波 动 较 大 的 值 ， 实 现 精准 定位 。 本 文 提出 数据 平滑 处 理 定义 : 
改进 平滑 算法 的 实现 分 为 两 个 阶段 :距离 的 估计 ;位 置 的 估计 。 Rn 是 第 n 个 平滑 估计 的 CSI 的 范围 ; 
2.2.1 距离 的 估计 Rw 是 第 个 预测 的 CST 的 范围 
本 文 提 到 的 CST 定位 机 制 由 参考 节点 和 未 知 节点 两 部 分 组 Ron 是 第 二 个 测量 的 的 CST 的 范围 ; 
成 的 。 在 二 维 坐标 系 中 表示 出 由 固定 位 置 的 参考 节点 的 位 置 坐 Raiw 是 第 个 平滑 估计 的 范围 比率 ; 
标 和 未 知 节点 收 到 的 CSI 值 进行 计算 。 求 二 者 的 标准 化 欧 氏 距 R prera 是 第 n 个 预测 范围 比率 ; 
离 。 4 ,6 表示 增益 常量 ,7 是 变化 的 时 间 段 。 
欧 氏 距离 计算 : 取出 参考 点 的 位 置 坐标 为 横 轴 的 值 记 做 并 s 处 理 步骤 如 下 : 
每 次 收集 到 未 知 节点 的 CSI 值 记 为 y; 则 三 者 的 标准 化 欧 氏 距 引 计 算出 Row) 和 Rao， 如 式 (17)(18) 所 示 。 
离 如 式 (15) 所 示 。 Resen =R pre AG, 7 Ryu) (17) 
d(x, y) = 4 — X)! * G5 7 y) (c, - y, Y Rs mV ER R pren) (18) 
b) 通 过 比较 第 nn 个 CSI 值 的 范围 与 比率 的 估算 和 第 n+1 个 
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CSI 范围 与 比率 的 估算 ， 如 式 (19)(20) 所 示 。 


c) 


R sretna) = Rao) * La (19) 
V rete = Voi (20) 
式 (20) 得 到 第 n+1 个 参考 点 的 值 ,再 由 式 (5) 计 算出 CSI 


的 值 ， 
以 CSI 
得 到 较 


2.22 位 置 的 估计 


传 


的 交点 ， 


确 精 确 


针 
文 提 出 


于 计算 过 程 中 加 上 了 传播 因子 对 信号 衰减 的 影响 ， 所 
值 存在 一 定 误差 , 把 求 得 的 CSI 值 进行 求 平均 数 的 处 理 ， 
为 准确 的 值 发 送 到 未 知 节点 。 


统 的 三 
如 图 2 
到 一 点 。 


边 测量 法 过 程 是 求 3 个 已 知 半 


所 示 ， 在 测 距 时 ，3 4p 


EA 
D 


图 2 ”三 边 定位 原理 医 
对 上 述 问 题 ， 为 了 减少 估计 坐标 与 实际 坐标 的 误差 ， 本 


£i 
I 
i" 
> 
HN 
I 
IL 
i. 
S 


标 圆 心 的 圆 
列 心 连接 线 无 法 准 


[rm 

IL 

cr 

[rii 

È 
| 


知 节 点 


a) M i p rp ax n ZB bz A s QS: y) 、 


了 改进 的 三 边 测 量 方法 ， 使 用 加 权 最 小 二 乘法 来 估计 未 
的 坐标 。 有 具体 方法 如 下 : 


(Gy), c7, 


(X44) ; 
b) 对 二 维 坐标 系 中 的 这 元 个 点 ， 用 处 于 样本 数据 的 中 心 位 


置 的 曲 


c) 引 入 平方 损失 函数 。 


线 来 拟 合 ; 


Q 为 残 差 平 方 和 ， 样 本 下 


HT 


FERAI n 


XQ DI. 


其 中 : 


则 通过 


变量 ， 


可 以 通 


Y,- B+ B, X, +e, Q1) 


e; =Y- ĝo- B, X, Q2) 


e 为 样本 Qo. y) 的 误差 。 平 方 损失 函数 为 


n 2 n X 2 n 2" m 2 
Q=} e, =), -Y:) =}, -2-2 X) Q3) 
i=l i=l i=l 


Q 最 小 确定 这 条 直线 ， 即 确定 入 po App A 


把 它们 看 做 是 2 的 函数 ， 就 变 成 了 一 个 求 极 值 的 问题 
过 求 导数 得 到 。 求 8 对 两 个 待 估 参 数 的 偏 导数 : 


99 ay X; - À,- P. XX-D-0 
op, f 

99 oy Y,- À,- B, X)-X) -0 
op, f 


Q4) 


函数 的 极 值 点 为 偏 导 为 0 的 点 。 即 活动 点 位 置 ， 解 得 : 


a _ nD XY- X -oY 
PTT A XP-AXY 

~ XY -XXY e 
Pi nYX/-Qxy 
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所 得 B. B, 即 为 所 求 活动 点 的 位 置 坐标 。 
2.3 E 
基于 上 述 理论 模型 分 析 ， 本 文 提出 基于 CSI 信号 的 人 员 轨 
m ic A 
开始 


skes 
[3 


———— "— —————— $— 


Music 算 法 估计 | 
FEDI 


三 边 定位 确定 未 
知 点 未 知 


3 轨迹 追踪 的 流程 图 

a) 通 过 布设 的 实验 环境 提取 CSI 信号 值 。 

b) 首 先 利用 Music 算法 估计 多 普 勒 频 移 ( 即 提出 的 Doppler- 
music 算法 )， 然 后 估计 AOA 频谱 ， 从 而 确定 人 员 位 置 是 否 在 
预 设 范围 内 ， 如 果 在 就 执行 步骤 c)， 否 则 转 回 步 又 a) 运行 。 

co) 采用 平滑 算法 对 步 台 b) 中 得 到 的 活动 点 进行 去 除 见 余 点 
里 ， 最 后 用 改进 的 三 边 定位 算法 确定 活动 点 位 置 ， 判 断 结果 
是 否 在 参考 范围 内 ， 如 果 在 就 执行 步骤 d)， 否 则 返回 。 

dj) 对 得 到 的 在 活动 范围 内 的 点 进行 处 理 ， 采 集 足 够 多 的 点 
模拟 出 轨迹 图 像 ， 程 序 结束 。 
2.4 通信 开销 和 算法 复杂 度 分 析 

为 了 计算 通信 开销 5 引入 一 下 几 个 参数 : N 表示 网 络 中 节 
i. A 表示 网 络 中 锚 节 点 的 数量 ，G 表示 网 络 的 平均 连 
表示 平均 邻居 节点 的 数量 ;及 表示 参与 一 定 多 边 测量 
定位 的 锚 节 点 个 数 。 

本 文 算法 的 开销 主要 表现 在 节点 间 的 通信 开销 和 数据 包 的 
传递 过 程 ,所 以 节点 间距 离 对 通信 开销 的 大 小 有 着 重要 的 影响 。 


E 
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当 未 知 节点 获得 与 3 个 或 3 个 以 上 锚 节 点 的 距离 , 则 执行 三 边 测 
量 定 位 。 设 未 知 节点 坐标 为 Ax, y), ， 锚 节点 坐标 为 五 (0 六) 7, 


H 


L, Gt. Y.) RATE A 8 88813 A E BAAN nhot DUI RI E 
立 线 性 方程 组 并 表示 为 式 (26) 的 形式 。 
Cx =b Q6) 


es = (yi —yx) 


其 中 : C=-2x Qn X) (i - Y LM e 
: y 
(KNX) Yay) 
m-n-Xx*x-x*x 
B- 5n-h xtX yl-yl T 


2 2 2 2 .2 2 
faf T Xka t Xy — Yk- t Yk 


建立 线性 方程 组 之 后 ,可 以 通过 最 小 二 乘法 求解 未 知 节 点 
的 位 置 估计 : 


x=(CTC) Ch Q9) 


首先 由 式 〈29) 得 到 未 知 节 点 的 位 置 ， 通 过 节点 位 置 计算 
相 邻 节点 之 间 的 距离 设 为 ,所 以 C 个 邻居 节点 的 平均 距离 为 


CX, Jie X =È Gx) 。 其 次 数据 包 的 传递 过 程 也 带 来 一 定 的 


通信 开销 ， 由 于 采用 可 控 的 泛 洪 在 网 络 中 传送 消息 ,每 个 锚 节 
点 发 送 广播 数据 包 , 中间 节 点 只 发 送 未 发 送 的 数据 包 , 每 个 节 
点 平均 发 送 A 个 数据 包 ， 因 为 要 经 过 两 次 泛 洪 过 程 ， 所 以 网 络 
中 每 个 节点 平均 发 送 2A 个 数据 包 ,N 个 网 络 节 点 的 通信 开销 为 
2AN。 故 总 的 通信 开销 E=(CX+2AN)。 

以 上 三 边 定 位 测量 得 到 未 知 节点 位 置 以 后 ， 使 用 最 小 二 
乘法 处 理 矩 阵 乘法 的 过 程 中 需要 锚 节 点 的 数量 K 和 参与 多 边 测 
量 定位 的 锚 节 点 个 数 A 分 别 计算 , 所 以 本 文 算法 的 时 间 复 杂 度 
为 T(AK)， 即 近似 为 T(n?)， 空 间 复杂 度 为 O Cn) 。 


3 ”仿真 实验 与 分 析 


3.1 实验 环境 
目前 ，CSI 信号 可 以 通过 Intel 5300 和 Atheros 9380 型 号 无 
线 网 卡 获取 ， 本 文采 用 Atheros 9380 网 卡 的 方案 ， 定 位 算法 所 
需 设 备 为 : 两 台 安装 有 Atheros 9380 网 卡 的 台式 电脑 ，CPU 型 
号 为 Intel Core i3-4150， 操 作 系 统 为 Ubuntu10.04LTS， 且 内 核 
和 无 线 网 卡 驱动 均 经 过 定制 ,其 中 一 台 台 式 机 作为 信号 发 送 机 ， 
另 一 台 作 为 接收 机 ， 实 验 场地 选择 9m X6m IZ AD, 
ET 25 个 正方 形 区 域 ， 每 个 正方 形 区 域 为 0.8mX0.8m， 接收 
机 与 发 射 机 相距 4.5 m， 天 线 高 度 为 1.2 m， 其 平面 图 及 详细 区 
域 划分 布局 如 图 4 所 示 ， 图 5 为 测试 人 员 现场 测试 图 。 
办 公 室 环境 的 平面 图 如 图 所 示 ， 实 际 测 试图 如 图 6 所 示 。 
本 文 的 跟踪 感应 区 的 大 小 是 9 米 X6 米 。 在 地 板 上 放置 一 些 标 
记 ， 并 使 用 摄像 机 记录 人 物 何 时 走 过 标记 。 在 实验 中 ， 选 择 5 
名 同学 在 2 周 的 时 间 内 (每 天 相同 的 时 刻 ) 沿 着 不 同 的 线条 ， 


GG 
o 


Cn 1 


矩形 ， 圆 等 轨迹 行走 。 本 文 为 每 个 人 收集 100 个 轨迹 ， 报 告 跟 
踩 错误 显示 追踪 的 性 全 
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图 6 办 公 室 环境 平面 图 
图 7 是 在 实际 环境 中 使 用 本 文 算法 实现 的 跟踪 结果 ; 


述 的 是 人 体 速度 方向 的 准确 性 。 


Oo 


em ; 
Q uet? El 0.8 
-10. 
-20. NE 
Lo fal 
-30 K 0.4} 
-40. f 
-50+ 0 
0 a ^ ——— Doppler-Music 
: 20 - . 一 一 一 'Kmeans 
50 7 0 00 20 40 60 80 100 
-20 Relative speed error 
图 7 轨迹 跟踪 结果 图 图 8 人 体 速度 方向 准确 性 
由 图 7 可 以 看 出 当 被 测 人 员 的 运动 轨迹 为 圆 形 时 ， 追 踪 效 


果 与 预期 结果 最 为 接近 ， 说 明 此 时 的 追踪 效果 最 好 。 由 图 8 说 
明 ， 中 速 运动 的 振幅 误差 小 至 14%， 而 方向 误差 仅 为 7”。 又 
将 本 文 所 提出 的 Droppler-Music 方法 与 基于 划分 的 K-means 算 
法 进行 比较 ， 由 于 K-means 算法 是 根据 CSI 的 幅度 来 估计 多 普 
勒 速度 ， 所 以 提供 不 了 方向 信息 ， 在 比较 过 程 中 ， 只 使 用 K- 


means 方法 来 估计 多 普 勒 速度 的 幅度 ， 并 使 用 Droppler-Music 
方法 来 提供 方向 信息 。 如 图 9 所 示 ， 在 加 上 方向 信息 来 提高 性 


能 的 前 提 下 ，K-means 方法 的 中 速 误差 仍然 为 44%， 方 向 误差 
为 17”， 远 大 于 Droppler-Music 方法 ， 原 因 是 Droppler-Music 


方法 估计 多 普 勒 速度 的 CSI 相位 ， 这 比 K-means 方法 中 使 用 的 
CSI 幅度 更 稳定 。 另 外 , 和 其 他 方法 相 比 ,本 文 的 Droppler-Music 
方法 可 以 同时 估计 多 普 勒 速度 的 幅度 和 方向 。 

天 线 功 率 (PA) 和 静态 分 量 (RS) 对 轨迹 追踪 精度 的 影响 如 
9 所 示 。 


aa 
图 9 相关 误差 图 

从 图 9 中 可 以 看 出 ， 当 不 调整 天 线 的 功率 也 不 去 除 静 态 分 
量 , 则 速度 误差 为 37%% ,方向 误差 为 26” ;只 调整 天 线 的 功率 ， 
相对 速度 误差 为 19%, 方向 误差 为 17”; 只 去 除 静 态 分 量 ， 相 
对 速度 误差 为 17%， 中 值 方向 误差 为 8”; 调整 天 线 功 率 又 去 
除 静 态 分 量 ， 相 对 速度 误差 为 11%， 中 值 方 向 误差 为 7”。 可 
见 ， 除 去 强 干 扰 静 态 分 量 后 ， 追 踪 性 能 显 着 提高 。 同 时 ， 调 整 
两 个 天 线 的 功率 ， 得 到 的 多 普 勒 速度 的 方向 信息 更 加 准确 。 
3.2 ”跟踪 性 能 分 析 

为 验证 分 组 丢失 /延迟 对 跟踪 性 能 的 影响 , 通过 控制 分 组 比 
较 丢 失 和 时 延 大 小 进行 实验 ， 结 果 如 图 10 所 示 : 
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图 10 
图 10 表明 ， 在 没有 丢 
误差 在 8% 左 右 , 在 实际 环境 中 , 由 于 环境 噪声 和 干扰 
的 影响 ,， 因 此， 随机 丢弃 一 些 CSI ži 
包 损 失 /延迟 情况 , 即使 有 50% 的 数据 包 丢失 , 本 方法 仍然 具有 
和 没有 数据 包 丢 失 相似 的 高 跟踪 精度 。 图 
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Xx (OD 


AURR 


昌 包 来 模拟 数据 


包 / 延 迟 的 情况 下 ， 该 算法 的 跟踪 


会 受到 


10 显示 了 在 在 实验 
室 中 跟踪 错误 的 CDF， 因 为 实验 室 遮 挡 物 较 多 ， 所 以 多 径 效应 


说 明 静 态 环 


比较 明显 ， 可 见 该 方法 也 能 实现 很 好 的 追踪 性 能 ， 
境 下 的 追踪 性 能 明显 高 于 动态 环境 下 。 

为 了 进一步 研究 速度 对 追踪 性 能 的 影响 ， 让 目标 以 三 种 不 
同 速度 步行 ， 慢 速 (<lm/s)， 中 速 (1~1.5Sm/s) 和 快速 (1.5~3my/s) 。 
速度 对 追踪 性 能 的 影响 如 图 11、12 所 示 。 


CDF 


10 


图 11 


速度 差异 影响 


11 说 明 目 标 的 速度 对 于 本 方法 的 跟踪 和 


20 
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速度 方向 误差 


图 12 速度 方向 影响 图 


FE 能 影响 ,3 


度 下 人 员 运 动 速 度 方向 对 跟踪 性 能 


13 所 示 : 


目标 间 能 达到 一 致 的 


慢 比 速度 最 快 的 时 候 误差 约 
好 ; 图 12 说 明 不 同 运 动 速 
的 影响 ， 说 明 速度 方向 误差 随 着 运动 速度 的 增 大 而 增 大 。 
本 方法 对 5 个 目标 进行 跟踪 ， 误 差分 析 如 图 
图 13 追 踩 误差 图 
由 图 13 可 以 看 出 ， 本 文 方法 在 不 同 
精度 。 
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踪 精 
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综 上 ， 通 过 实验 分 析 得 到 ， 在 丢 包 率 、 速 度 差 异 、 方 向 差 


多 人 效应 都 存在 的 情况 下 ， 本 文 算法 在 定位 方面 


。 有 效 的 缓解 
度 的 问题 。 


了 


有 很 大 的 


Sn 


于 多 径 衰弱 、 非 视 距 环境 影响 带 来 的 追 


好 的 解决 思路 。 


尤其 是 在 速度 方向 误差 方面 本 文 算法 提出 了 一 


录用 稿 


3.3 ”性 能 比较 

为 了 进一步 表现 本 文 算法 在 追踪 性 能 方面 的 优势 ， 与 传统 
经 典 算法 K-means 和 PCA 算法 进行 对 比 。 通 信 开 销 性 能 对 比 
14 所 示 ， 算 法 复杂 度 性 能 对 比如 图 15 所 示 。 
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图 14 算法 通信 开销 对 比 图 
图 14 可 知 ， 本 文 算法 随 着 时 间 的 推移 更 早 的 表现 出 较 
稳定 通信 开销 ， 且 明显 优 于 PCA 算法 。 性 能 越 稳 定 表 示 通 信 开 
销 越 低 ， 得 到 的 节点 位 置 就 越 准 确 ， 所 得 轨迹 追踪 的 性 能 就 越 
好 。 由 图 15 可 知 ，Doppler-music 算法 的 时 间 复 杂 度 在 一 段 时 
间 后 始终 保持 不 变 且 性 能 明显 高 于 其 他 两 种 算法 , 说 明 较 PCA 
P a A N RE E ER 
销 的 同时 降低 能 耗 ， 从 而 实现 更 加 精准 的 追踪 轨迹 。 


2 3 
算法 时 间 复 杂 度 


图 15 算法 时 间 复 杂 度 对 比 
综 上 ， 基 于 Doopler-Music 算法 的 CSI. 轨迹 追踪 算法 在 通 


言 开 销 和 时 间 复 杂 度 上 都 明显 优 于 其 他 两 种 算法 ， 为 实现 精准 


轨迹 奶 踪 提出 一 种 新 的 解决 思路 。 


4 tmi 


本 文中 ， 旨 在 设计 一 个 利用 CSI 信号 实现 室内 人 员 轨 迹 追 
踪 的 方法 。 实 现 了 一 个 人 的 轨迹 精确 追踪 。 以 后 的 工作 中 需要 
解决 的 问题 有 多 人 追踪 、 实 时 追踪 等 ， 室 内 人 员 方 法 的 应 | 
有 很 广泛 的 前 景 。 此 外 一 个 重要 难点 还 有 目前 的 实验 环境 是 
2.4GHZ 的 ， 下 一 步 工作 就 是 把 开发 环境 移植 到 SGHZ 的 环境 
下 ， 以 获得 更 加 精准 的 定位 性 能 。 
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